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Abstract 

Using perturbation theory in the strong coupling regime, that is, the dual Dyson series, and 
renormalization group techniques to re-sum secular terms, we obtain the perturbation series of 
the two- level system driven by a sinusoidal field till second order. The third order correction to 
the energy levels is obtained proving how this correction does not modify at all the localization 
condition for a strong field as arising from the zeros of the zero-th Bessel function of integer 
order. A comparison with weak coupling perturbation theory is done showing how the latter is 
contained in the strong coupling expansion in the proper limits. The strong coupling expansion 
we obtain proves to be accurate in the regime of high-frequency driving field. This computation 
gives an explicit analytical form to Floquet eigenstates and quasi-energies for this problem, for 
high-frequency driving fields, supporting recent theoretical and experimental findings for quantum 
devices expected to give a representation for qubits in quantum computation. 

PACS numbers: 72.20.Ht,73.40.Gk,03.67.Lx 
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I. INTRODUCTION 



Quantum computatio n 1 ' 2 ^ 4 , intended as information theory performed using qubits, that 
is, two-level systems evolving by the unitary evolution as obtained by the laws of quantum 
mechanics, demands an increasing control on such systems. 

Recent experimental findings^^^^ proved that qubits can be realized by solid state 
devices, making the realization of a quantum computer even more possible. Such solid state 
devices performs in a way or another like two-level systems driven externally by some time- 
varying field. This means that it is very important a clear theoretical understanding of them 
to make their control easier. 

Advances in this way have been realized by devising new approaches to the solution of 
the Schrodinger equation in different regime a 1( | )i1111211 ? 11411 ^ . These methods are perturbative 
in their nature but permit to study a quantum system in different physical regimes. So, 
they seems to be a proper framework to make two-level systems theoretically manageable 
in any situation. These methods give to the Floquet method, for periodical perturbations, 
a strong analytical support. 

The proper control of a qubit, realized by some solid state device, relies on the possibility 
to obtain localization of a particle in each one of the two available states. This effect is 
also known as coherent destruction of tunneling (CDT) 16 ' 17 as, tunneling between one state 
and the other is destroyed by making unity the probability of staying of the particle in 
the initial state. So, this question has been recently the main motivation for several works 
on two-level system o 18 ' 19 i 20 i 21 . Some of these papers have put forward an interesting non 
perturbative result showing a set of curves for the crossing of the quasi-energies of the 
Floquet solution for two-level systems with different periodical external fields. As it stands, 
such a result, being not perturbative in its very nature, appears rather difficult to recover 
with perturbation theories. 

The aim of this paper is to show how perturbation theory can give definite results for 
localization in different regimes, putting the study of two-level systems on a sound ground 
from an analytical point of view. Besides, an expression is given for localization, after 
that an effective time-independent Hamiltonian is obtained using perturbation theory. Our 
results confirm till third order the well-known result that the zeros of the zero-th order 
Bessel function are the points where CDT occurs. Besides, our result recovers the small 
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perturbation theory computation in the proper limits. This is in agreement with the recent 
numerical computations^^. But, it is essential to emphasize that the main result of the 
paper, the strong coupling expansion, is very accurate for high-frequency driving fields. 

The paper is so structured. In Sec. |H]we discuss localization on a quite general ground, 
using Floquet theory. In Sec. 11111 we give a presentation of the perturbation methods we 
use to obtain our results. In Sec. |lV|we present the computation for the two-level system 
driven by a sinusoidal field giving the effective time-independent Hamiltonian obtained by a 
computation till second order. In Sec. Elwe give the third order correction to the effective 
Hamiltonian showing how the localization condition at high fields is kept in agreement with 
all numerical results as also happens to the small field condition. Then we discuss the 
conditions under which localization happens comparing the perturbation series obtained 
with the weak perturbation series showing the way one result contains the other. Finally, 
in Sec. IVII conclusions are given. 



II. LOCALIZATION IN DRIVEN TWO-LEVEL SYSTEMS 

The model we consider has the Hamiltonian 

H=^-a 3 + ga 1 f(t) (1) 

being A the separation between the two levels, g the coupling, o\ and a 3 Pauli matrices 
and f(t) = f(t + T) with T the period of the perturbation. The theory of CD T 16 i 17 can be 
straightforwardly applied. Given the unitary evolution operator U(t), one can compute the 
probability of the system of being in the initial state as 

p(t) = \(m\u(t)\m)\ 2 (2) 

with |^(0)) being one of the eigenstates of a 3 . One has CDT when this probability is equal 
to unity. As the Hamiltonian is periodic in time, one can apply the Floquet theory and 
take22 

P = -[ P(t)dt (3) 



T jo 

to ascertain that we have true localization by having unity also in this case. The consequence 
is that one can show, again from Floquet theory, that the conditions for CDT arise from 
crossing of quasi-energies. This condition is necessary but not sufficient. 



Our aim in the following is to derive such conditions that apply to the quasi-energies for 
the sinusoidal case, that is, assuming f(t) = cos(ut) with to = 2n /T. We give an explicit 
analytical expression for the Floquet eigenstates and quasi-energies by using perturbation 
theory in the strong coupling regime with the method described in the next section. 

III. DUAL DYSON SERIES AND RESUMMATION TECHNIQUES 

It is generally believed that perturbation theory, being plagued by secularities , that is, 
unbounded terms in a perturbation series that increase as powers of time, t,t 2 , ■ ■ •, should 
be avoided to treat a model like the one we consider in this paper. Indeed, the meth- 
ods devised so far for removing such singular terms proved to be generally not very easy 
to apply. Recently, a new approach, that can be called dynamical renormalization group 
metho d 23 i 24 i 25 i 26 i 27 i 28 , has given an algorithmic way to remove secularities, making computa- 
tion in perturbation theory just tedious but very easy to accomplish. This method, coupled 
with the dual Dyson serie o 29l3 ° l31l32lS3 , gives a straightforward method to compute higher or- 
der corrections to well known results^. Besides, we will get an explicit analytical expression 
for the Floquet modes and the quasi-energies in the given approximation and we will prove 
that the small coupling result is recovered. This is a property of the dual Dyson series. 

The dynamical renormalization group method was firstly formulated by the Urbana 
Group 23 ' 24 based on the observation that renormalization group methods can be seen as 
a means of asymptotic analysis^ 3 -. They showed how their approach can be proven equiva- 
lent to other methods that resum secularities, such as the multiple time scale method 3 - 4 -, for 
many cases. This method, for its very nature, can be described by an example. So, let us 
consider the well-known equation of the forced harmonic oscillator 

x(t) = -x(t) - ex{tf (4) 

being e the strength of the anharmonic term assumed to be small. A naive perturbation 
expansion in e till first order gives the well-known result 

3 A 
x{t) = A sin(i + O ) + e-A (t - t ) cos(t + O ) - sin(3t + 30 o ) + 0(e 2 ) (5) 

8 lb 

being A and 0o constants depending on the initial conditions. We recognize a secular term 
t — to that makes this series useless as it breaks down for e(t — to) > 1- The situation 
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can be improved if we interpret the time to as the logarithm of the ultraviolet cutoff in 
quantum field theory and introduce A and as the renormalized counterpart of Aq and 
0o due to the fact that the nonlinearity may change these constants. Then, we introduce 
another renormalization point by splitting t — t in t — t + r — t and adsorb the terms 
containing t — t into A and 0. We introduce a multiplicative renormalization constant 
Z\ = 1 + Y^Li a n£ n and an additive one Z 2 = Y^Li b n £ n so that A Q (to) = Zi(t , t)A{t) and 
0o(^o) = 0( r ) + Z2{t ,r) with the coefficients a n and b n to be computed order by order to 
remove the terms with t — to as happens in standard renormalization group 35 i 36 . It is easily 
seen that, in our case, a possible choice to first order is ai = and b± = — §(t — to) removing 
the secular term and we are left with 

3 A 
x(t) = A sin(t + 0) + e-A(t - r) cos(t + 0) - e— sin(3t + 30) + 0(e 2 ). (6) 

8 16 

being now A and function of r. But r does not appear in the original problem and so x(t) 
must be independent on it. This gives us the condition 



dx 



(7) 



T=t 



that is the renormalization group equation that in our case gives the equations 

where the well-known shift in the frequency of the oscillator is recovered. Besides, the 
secularity is completely removed by taking into eq.© the condition r = t with A(t) and 
0(t) given by the solutions of the renormalization group equations. We have finally 



x(t) = A(Q) sin 



l + ~eA{0))t+ <f)(0) 



A(0) . 



3(l + -eA(O))t + 30(O) 



+ 0(6 2 ). 

(9) 

So, a straightforward application of renormalization group methods has permitted to remove 
a secular term in the perturbation series. But there is a way to make the computation simpler 
by noting that we have done nothing else than to compute the envelope of eq.© using known 
renormalization techniques. 

Then, the method of renormalization group to resum secularities in a perturbation series 
that we present here is obtained by the mathematical theory of envelopes and is due to 



Kunihir o 25 i 26 i 27 . In order to describe it, let us consider the following equation 



x{t) = f{x{t),t), 



(10) 



with x(t) that can be also a vector. The initial condition is given by x(to) = X(t ). At 
this stage we assume X(t ) not yet specified. We write the solution of this equation as 
x(t; to, X(to)) that is exact. If we change to to t' we are able to determine X(to) by assuming 
that the solution should not change 



x(t;t ,X(t ))=x(t;t ,X(Q) 



that in the limit to ~~ > t' Q becomes 



dx dx dx dX 
dt dt dX dt 







(11) 



(12) 



giving the evolution equation or flow equation of the initial value X(to). We see again the 
renormalization group equation proper to this approach. 

Till now, all our equations are exact and no perturbation theory entered in any part 
of our argument. But, except for a few cases, the solution x(t; to, X(t Q )) is only known 
perturbatively and such a solution is generally valid only locally, i.e. for t ~ t and t ~ t' 
and a more restrictive request should be demanded to our renormalization group equation 

dx dX 



dx 
dto 



t =t 



dx 



+ 



to=t 



dX Otn 



0. 



(13) 



t =t 



But this equation can be interpreted by the mathematical theory of envelopes 2 ^. Indeed, 
varying to we have that x(t; t , X(t )) is a family of curves with t as a characterizing 
parameter. Then, eq.(|13|) becomes an equation to compute the envelope of such a family 
of curves. Such an envelope is given by x(t; to = t) = X(t), the initial condition. It can be 
proven that X(t) satisfies the equation (fTU|) in a global domain up to the order x(t;to) for 
t ~ t . 

Now, we can describe our approach in some steps to show how such an algorithmic 
perturbation method indeed works 

1. Consider the following unitary transformation on the Hamiltonian Q (here and in the 
following we set % — 1) to remove the perturbatio n 29 ' 30 

rt 



Up(t) = exp 



-'.(/"" l / f(t')dt' 



(14) 



giving the transformed Hamiltonian 

h f(P) = y^exp 



-i2ga x f f(t')dt' 
Jo 



the dual Dyson series is computed p y 29 i 30 i 31 i 32 

SD(t,t ) = Texp 



-%e I H F (t')dt' 

to 



(15) 



(16) 



being as usual T the time ordering operator and an ordering parameter e has been 
introduced that will be taken unity at the end of computation. It is fundamental for 
our argument that the computation of this series is performed at a different starting 
point to- 



2. Assume, at the start, that the time evolution operator has the form 

U(t,t ) = U F (t)S D (t,t )U R (t ) 
where U R (t ) is a "renormalizable" part of the unitary evolution. 

3. At the given order one gets S£>(t, to) as 

S D (t,t Q ) =I-iefi(t,to) -e 2 / 2 (Mo) + ••• 



(17) 



(18) 



and, at this stage, if some oscillating functions in t appear like e~ luJt ° then introduce 
the phase 0(to) = —to as a "renormalizable" parameter rewriting it as e""^''^ We 
make this choice assuming that any initial phase of the system can be changed by the 
dynamics. The minus sign is fixed arbitrarily. The secularities must be left untouched. 



4. Eliminate the dependence on t by requiring 25 

dU(t,t ) 



dtr 







t =t 



and one obtains the renormalization group equation 

dU R (t) 



dt 

d<p{t) 
dt 



egiU R (t) + e 2 g 2 U R (t) + 0(e 3 ) 
e0 1 0(t) + e 2 2 0(t) + O(e 3 ) 



(19) 



(20) 



where, at some stage, to obtain such equations at the second order, we have to use 
their expressions at the first order and, to compute their form at order n-th, one has 
to use these equations at the order {n — l)-th, into the condition (fT9*|) . This is a step 
toward the computation of the envelope of the perturbation serie o 25 ' 26 . 



5. Finally, the renormalization equations should be solved and substituted into the equa- 
tion 

U(t,t )\ t0=t (21) 

giving the solution, i.e. the envelope, we were looking for without secularities at the 
order we made the computation. 

We will give an explicit application of this procedure in the next section choosing f(t) = 
cos(u;t), that is, a sinusoidal driving. 

IV. STRUCTURE OF THE PERTURBATION SERIES TO SECOND ORDER 

The unitary transformation {7p(t) for a sinusoidal driving field takes the form 



U F (t) = exp 



— i— sin(o;t)(7i 



and the transformed Hamiltonian takes the form 



A 



H F (t) = — c^exp [— iz sm(ivt)ai] . 



(22) 



(23) 



having put z — -A This Hamiltonian gives the term <Sd(£, to) till first order as 

S D (t,t ) = I - i^a 3 J (z)(t - t ) (24) 

/\ g— inwta\ g— inuitgcri 

+ 17^1 J n( Z ) h ' • ' 

2 ~f n mo 

where J n (z) are the Bessel functions of n-th order being n an integer. We see immediately 
that a secular term, proportional to t — 1 , plagues our computation. To remove it we rewrite 
the above expression as 



S D (t,t ) = I — i—a 3 J (z)(t — t ) 



(25) 



+—cr 3 cri J n{z) V ... 

2 nUJ 



introducing into the last oscillating term the renormalizable parameter 0(to). So, one has 



dU(t,t ) 
dtr, 



U F (t) 



+ U F (t) 



.A 



a 3 J (z) - —a 3 Jn(z)e ir 



oMtoW d <j>(to) 



dto 



U R (t ) 



1 ~ i-^0zJo(z)(P ~ *o) + ^-<T 3 <Ti Y J n{ z ) 1" • • ■ 

z i n _^ nuj 



(26) 

dU R (t ) 
dto 
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and imposing the condition (|19j) we get the renormalization group equations 



dU R (t) 
dt 

d(j)(t) 
dt 



-i^Mz)U R {t) + O 



A x 2 



(27) 



O 



A x 2 



where use has been made of the fact that these equations are at least first order, transforming 
some terms in eq. ()26|) into second order, making them negligible. We see that we have 
recovered a well known result of the series for high frequency^ with its corrections to first 
order. This gives the final result, using eq. pTf . 



U(t,0) = U F {t) 



I + 77^3^1 V Jn{z) h 

2 ^ no; 



-iyo-3 Jo(z)t 



(28) 



This result is well-known in literatur o n i 12 i 31 i 33 . It is also been proven that this is the expres- 
sion of the Floquet unitary evolution at this order with the given approximation ^ 1 1 > 12 . We 
see that we have a periodic part and a part that originates the quasi-energies. 

So, what is really interesting is to get higher order corrections. Particularly, we see that 
the quasi-energies are given by ±^Jo(z), a well-known fact, so that, CDT occurs, in this 
approximation, at the zeros of the Bessel function Jq(z)H. We now prove that higher order 
corrections preserve such a result. 

Applying the above procedure till second order it is straightforward to obtain 



U(t,0) = U F (t) 
A 2 



A 



-incutai 



l 



nuo 



(29) 



i—oxMz) ^ J n {z) 
2 ^ n 2 u 2 



A 2 



A 2 



Arwjt 



n 2 u 2 



E 



J n\ {.Z)Jri2 \Z) 



n 2 (ni - n 2 )u; 2 



711712^ 



if <7 3 Jo(z)t+i| ;r F(z)o-iJo(z)i 



where we have introduced the function F(z) = Y, n ^o the last exponential. We 

recognize a product of a periodic unitary operator and a term originating quasi-energies. 
So, it is interesting to see that now we have an effective Hamiltonian giving the quasi- 
energies with a correction term proportional to o\. But, the most important point is that 
this correction is again proportional to Jq(z) and so, the zeros of this Bessel function gives 
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CDT also at second order. The effective Hamiltonian can be written as 

= - (30) 

that can be easily diagonalized but we do not pursue this matter here. Indeed, our aim in 
the next section will be to give the third order correction to this Hamiltonian proving that 
the high-frequency CDT is again determined by the zeros of Jq(z). 



V. THIRD-ORDER CORRECTION TO THE QUASI-ENERGIES AND LOCAL- 
IZATION 



The algebra is rather tedious but straightforward and the application of our method gives 
the third order correction to the effective Hamiltonian 



H. 



(3) 
eff 



A 3 J (z) _ Jl(z) 



Au 2 



E 



n- 



-0-3. 



(31) 



This is a correction to the 03 term of the effective Hamiltonian so, we can conjecture that 
even order corrections go to the U\ term and the odd order corrections go to the cr 3 term. 
We see again that this term is proportional to Jo (z) confirming the exactness in the high- 
frequency limit of the occurring of the CDT at the zeros of such a Bessel function. Such 
correcting terms into the effective Hamiltonian can be recognized as a.c. Stark shifts and 
Bloch-Siegert shifts and are relevant, eventually, at the resonanc e 2 ^? where Rabi flopping 
also in this regime is expected with the renormalized levels ±^-Jo(z). 

Anyhow, a further check can be obtained with the small coupling perturbation theory. 
Our method can be applied again but now we use the interaction picture. So, if we introduce 
the unitary transformation Uj(t) = e~ l ^°" 3 *, we get the transformed Hamiltonian 



Hj = g cos(ujt)aie 



-iAa 3 t 



(32) 



and we can built the Dyson series, out of resonance (A 7^ uj), to second order as 



U(t,0) = U^t) 



g / e t(w-A)<7 3 t _ y 
2 V uj — A 



uj + A 



(33) 



g 2 /cos(2cjt) + ia 3 ^ sm(2ut) e i(w+A)<T3 * - 1 
+ T I uj 2 - A 2 to 2 - A 2 



+ 



{uj - A) 2 



-i(u}-A)a 3 t e i(u+A)a 3 t ,,2 

+ — -2 



^ + A 2 



A 2 



(u + Ay 



A 2 ) 2 



+ 
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A rather interesting result is that this series holds for any ratio — differently from the dual 
Dyson series discussed above that holds for - < 1. From the effective Hamiltonian (|30|) 
and the third order correction we see that the coefficient of the 03 part is given by 

AJ (z) A*J (z) J 2 (z) 
S3 -~2 (34) 

2 

that in the small coupling approximation, being J ( z ) ~ ¥ 1 — ir> an d considering just the 

2 

first term of the series giving 4j-> we obtain 

A z 2 A 3 z 2 

From the small coupling expansion we recognize, after a simple rearranging of the terms, 
the CT3 term 

A g 2 A 

s * = 2-^u=& (36) 

that, in the expected limit — -C 1, is the same as S3. Similarly, one can verify the same result 
at any order of the two perturbation series. This gives an a posteriori verification of our 
computations. Besides, it confirms our point that the dual Dyson series is a high-frequency 
series while the Dyson series holds for both the limits of high and low frequencies, but the 
result of the Dyson series in the high-frequency limit is contained in the dual Dyson series 
when the limit of small coupling z — > is also taken. Indeed, the Dyson series and the dual 
Dyson series can be taken to coincide in the limits of small z and - < 1, giving an interesting 
relationship between the transformed Hamiltonians we use to obtain the perturbation series. 

We are now in a position to discuss the CDT in the limit of small coupling z — > 0. For 
our aim, we consider the analytical expression of the quasi-energies obtained in Ref.^. This 
gives 

A\ 2 2 Asin(7r£)' 



e± = ± — arccos 
2tt 



COS 7T — + 7T 



, mod(o;) (37) 



ujj '"4a; 1 - ^ 

that yields, in the limit z < 1, eq.([36|) proving again the correctness of our computations. 
The important point is that eq.(|37|) contain the result for CDT that the quasi-energies 
crosses for — = 2n being n an integer and this cannot happen for the perturbative result 
unless we are able to resum the series. So, we can draw the conclusion that perturbation 
theory by Dyson series can prove to be very effective for the study of CDT for strong fields 
and high-frequency regime. 



II 



VI. DISCUSSION AND CONCLUSIONS 



From the von Neumann- Wigner theorem^ one has to expect that the crossings of the 
energy levels form a one- dimensional manifold. This result is central to fully understand 
CDT for the two-level system we considered. Recent works on this question derived such a 
result numericall y 20121 . In this paper we have shown that, unless some smart resummation 
technique is applied to a perturbation series, perturbation theory is helpful just to study the 
behavior of the model in different regimes, recovering CDT under certain conditions. 

The reason for this conclusion lies on the fact that the content of the von Neumann- 
Wigner theorem is non-perturbative. So, if we are able to resum a perturbation series to 
recover a non-perturbative result, there may be a possibility to give an analytical proof of 
the numerical results obtained so far. 

Notwithstanding such a conclusion, we have proved the power of the Dyson series in the 
study of two-level systems in regimes where other means may prove unsuccessful, giving 
explicitly the form of the Floquet modes till second order and quasi-energies till third order, 
for a strong field and high frequency. At this stage, it appears mandatory to introduce 
resummation techniques in perturbation theory to complete the algorithmic approach we 
discussed in this paper. 

Different regimes of parameter space have been analyzed by our perturbation technique. 
Notably, we have obtained perturbation series for the high-frequency regime A <C u in the 
strong coupling approximation g ^> A by dual Dyson series. The Dyson series has given 
the small coupling expansion g < A that holds for any ratio ^ but when this ratio is taken 
to be small, the dual Dyson series recovers the Dyson series for g -C A as we have shown. 
These results appear very promising for the study of quantum systems in different regimes 
by a general way to approach computations with perturbation theory. 

Acknowledgments 

I would like to thank Charles Creffield for suggesting me such an interesting question. 



* e-mail: marcofrasca@mclink.it 



12 



1 D. Bouwmeester, A. Ekert and A. Zeilinger (Eds.), The Physics of Quantum Information, 
(Springer, Berlin, 2000). 

2 C. H. Bennett and D. P. DiVincenzo, Nature 404, 247 (2000). 

3 R. Blatt, Nature 412, 773 (2001). 

4 D. Kielpinski, C. Monroe and D. Wineland, Nature 417, 709 (2002). 

5 T. H. Oosterkamp, T. Fujisawa, W. G. van der Wiel, K. Ishibashi, R. V. Hijman, S. Tarucha 
and L. P. Kouwenhoven, Nature 395, 873 (1998). 

6 Y. Makhlin, G. Schon and A. Shnirman, Nature 398, 305 (1999). 

7 Y. Nakamura, Yu. A. Pashkin and J. S. Tsai, Nature 398, 786 (1999). 

8 B. E. Cole, J. B. Williams, B. T. King, M. S. Sherwin and C. R. Stanley, Nature 410, 60 (2001). 

9 Yu. A. Pashkin, T. Yamamoto, O. Astafiev, Y. Nakamura, D. V. Averin and J. S. Tsai, Nature 
421, 823 (2003). 

10 J. C. A. Barata, and W. F. Wreszinski, Phys. Rev. Lett. 84, 2112 (2000). 

11 J. C. A. Barata, Ann. Henri Poincare 2, 963 (2001). 

12 J. C. A. Barata, D. A. Cortez, J. Math. Phys. 44, 1937 (2003). 

13 V. Delgado, and J. M. Gomez Llorente, J. Phys. B 33, 5403 (2000). 

14 A. Santana, J. M. Gomez Llorente, and V. Delgado, J. Phys. B 34, 2371 (2001). 

15 K. Fujii, |quant-ph/0301145| 



16 F. Grossmann, T. Dittrich, P. Jung and P. Hanggi, Phys. Rev. Lett. 67, 516 (1991). 

17 M. Grifoni and P. Hanggi, Phys. Rep. 304, 229 (1998). 

18 F. Grossmann and P. Hanggi, Europhys. Lett. 18, 571 (1992). 

19 C. E. Creffield and G. Platero, Phys. Rev. B 66, 235303 (2002). 

20 J. M. Villas-Boas, Wei Zhang, S. E. Ulloa, P. H. Rivera and N. Studart, Phys. Rev. B 66, 
085325 (2002). 

21 C. E. Creffield, Phys. Rev. B 67, 165301 (2003). 

22 J. Shao and P. Hanggi, Phys. Rev. A 56, R4397 (1997). 

23 L. Y. Chen, N. Goldenfeld and Y. Oono, Phys. Rev. Lett. 73, 1311 (1994). 

24 L. Y. Chen, N. Goldenfeld and Y. Oono, Phys. Rev. E 54, 376 (1996). 

25 T. Kunihiro, Prog. Theor. Phys. 94, 503 (1995); (Erratum), Prog. Theor. Phys. 95, 835 (1996). 

26 T. Kunihiro, Phys. Rev. D 57, R2035 (1998). 

27 S.-I. Ei, K. Fujii and T. Kunihiro, Ann. Phys. (NY.) 280, 236 (2000). 

13 



M. Prasca, Phys. Rev. A 58, 771 (1998). 
M. Frasca, Phys. Rev. A 45, 43 (1992). 
M. Frasca, Phys. Rev. A 47, 2374 (1993). 
M. Frasca, Phys. Rev. A 56, 1548 (1997). 
M. Frasca, Phys. Rev. A 58, 3439 (1998). 
M. Frasca, Phys. Rev. A 60, 573 (1999). 

J. Kevorkian and J. D. Cole, Perturbation Methods in Applied Mathematics, (Springer, New 
York, 1981). 

E. C. G. Stiickelberg and A. Petermann, Helv. Phys. Acta 26, 499 (1953). 

M. Gell-Mann and F. E. Low, Phys. Rev. 95, 1300 (1954). 

M. Frasca, J. Opt. B: Quant. Semicl. Opt. 3, S15 (2001). 

X.-G. Zhao, Phys. Rev. B 49, 16753 (1994). 

J. von Neumann and E. Wigner, Phys. Z. 30, 467 (1929). 



14 



